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Correlation energy of two electrons in a ball 
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We study the ground-state correlation energy of two electrons of opposite spin confined within a D- 
dimensional ball [D > 2) of radius R. In the high-density regime, we report accurate results for the exact 
and restricted Hartree-Fock energy, using a Hylleraas-type expansion for the former and a simple polynomial 
basis set for the latter. By investigating the exact limiting correlation energy Ec^'' = Mmn^o Ec for various 
values of D, we test our recent conjecture [J. Chem. Phys. 131 (2009) 241101] that, in the large-D limit. 
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for any spherically-symmetric confining external potential, where S = ^/iD — 1). 



PACS numbers: 31.15.ac, 31.15.ve, 31.15.xp, 31.15.xp, 31.15.xr, 31.15.xt 

Keywords: correlation energy, two-electron systems, spherium, hookium, Hylleraas expansion, Hartree-Fock 
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I. INTRODUCTION 

In the early days of quantum chemistry, there was con- 
siderable interest in cavity-confined atoms as a model 
for high-density atomic gaJ^^ and extrapolation of 
high-density results provides a convenient but powerful 
route to understanding the intermediate-density regime.^ 
Thanks to Hylleraas' worlP, the compressed helium-like 
ions have been widely studie(j^ and interest in these 
continues unabatedP^ Other confined systems such as 
electrons in squar^l^EH^ cylindrical^ and sphericaP^^t^ 
boxes have also attracted attention. The last of these 
has been extensively used for the assessment of density- 
functional approximations^^ and the study of Wigner 
molecules^^ at low densities. ^^ -^^ 

In a previous article,'22lwe studied the high-density cor- 
relation energy Ei°^ for various two-electron systems con- 
fined to a I?-dimensional space {D > 2) by an external 
potential V{r) oc r"\ As the high-density limit sheds 
light on intermediate densities, the large-dimen sion l imit 
provides useful insights into the D = 3 case.f^S^H Por 
the helium-like ions (m = —1), the spherium atomd^&HMI 
(to = 0), and the Hooke's law atom^^SHal (to = 2), we 
found that, in the large- 1? limit. 



(1) 



where S — 1/{D — 1) and the coefficient C w 1/6 varies 
slowly with to. On this basis, we conjectured that Eq. ([T]) 
is true for any spherically-symmetric confining external 
potential. 

At the end of our previous workj^^^l we observed that 
it would be highly desirable to consider D-ballium, the 
system in which the two electrons are trapped in a D- 
dimensional ball of radius R. This model is a severe test 
of our conjecture because it corresponds to m = oo. 
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The present study focuses mainly on the high-density 
regime {R > but small) and the corresponding lim- 
iting case (i? = 0). We report accurate results for the 
restricted Hartree-Fock (HF) and exact energies (Se c. [Til 
and III respectively). For the limiting case (Sec. |IV[ ), 
perturbation theory is used to expand both the HF and 
exact energies and this allows us to determine the limit- 
ing correlation energy in Z?-ballium. We use atomic units 
throughout. 

The Hamiltonian of Z3-ballium is 



H 



2 



1 



+ V{r^) + Vir2)+ , 
2 ri2 



(2) 



where ri2 = |ri — is the interelectronic distance, and 
the external potential is defined by 



V{r) 



0, if r < R, 
oo, otherwise. 



(3) 



Any physically acceptable eigenfunction of ^ must sat- 
isfy the Dirichlct boundary condition 



^(n = i?) = *(r2 = i?) = 0. 



(4) 



II. RESTRICTED HARTREE-FOCK APPROXIMATION 

The spin-restricted HF solution,'^ which is the only 
HF solution in the high-density regime, is given by 



(5) 



If we introduce the scaled coordinate t ~ r/R, the HF 
orbital (j){t) is an eigenfunction of the Fock operator 



F : 



2^^? + ^^^(^)- 



(6) 
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where 
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FIG. 1. Accuracy A of the HF energy of 3-ballium [R = 1) 
with respect to basis set size A''. TA results taken from Ref. 
[21] and LG from the present study. 



For S states in a Z?-diniensional spacej^^'^ we have 

D-l d 



df^ 



Iq max(r, xj 



where a = ™™(*''^) and F is the hypergeometric function. 

Unhke Thompson and Alavif^U who expanded the HF 
orbital in a basis of spherical Bessel functions,''^ we chose 
to explore an even-degree polynomial basis, writing 

N-l 

m = (1 - i') E ^fc*''- (9) 

fc=0 

Any such orbital is smooth at the center of the ball, i.e. 
(t>'{0) = 0, 

and is cusped and vanishes at the boundary, i.e. 

0'(1)<O, =0. 
It can be shown that the resulting HF energy is 

IT 1 U 



t dt' 
3-D 1 D 



(7) 
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x'^-^dx, 



(8) 



i?HF — 



S RS^' 



with 



CiCj 



CiCj 



D 



Aij 



(«+J + f)2 (*+J + f-l)3 
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(10) 
(11) 
(12) 

(13) 
(14) 

(15) 



(a)f, 



r(a + b) 

r(a) 



(16) 



is the Pochhammer symbol and F is the Gamma 
function.ES The coefficients /3„ and 7„ are given by 
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7ri 



iF2{ 
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l + f;l) 



(17) 
(18) 



where pFq is the generalized hypergeometric function.^!] 
For = 3, the hypergeometric functions reduce to unity. 

The energy ( 12 1 can be minimized with respect to the 
coefRcients Cfe using a numerical solver,!^ thus avoiding 
the self-consistent field procedure usually needed for this 
kind of calculation.l3SI 

Henceforth, we define the accuracy of an energy E as 



^ = -logio [{E-E*)/E*] 



(19) 



where E* is our best estimate of the exact energy. In 
loose terms, A is the number of correct decimal digits. 

Figure [T] shows how the accuracy of the HF energy of 
3-ballium (7? = 1) improves as N increases. For very 
small N , the spherical Bessel basiP^^ (TA) is more accu- 
rate than the polynomial (LG) basis. However, although 
both the TA and LG bases seem to exhibit exponential 
convergence as N increases, the TA energy improves by 
roughly one order of magnitude and the LG energy by 
roughly two orders of magnitude as each basis function 
is added. As a result, one obtains the HF energy to 20 
digits using ([9| with N ~ 10. However, the origin of the 
superiority of the polynomial basis is not clear. We find 
that the resulting expansion coefficients Cfc decay roughly 
exponentially and the convergence behavior for other D 
is similar. 

Numerical results for 3-ballium (i? = 1) are shown in 
Table |l] For N = 7 basis functions, the Bessel and poly- 
nomial bases yield HF energies that lie 13 ni^h and 2.8 
pi?h above the HF limit, respectively. Analogous behav- 
ior is observed for the larger values of R, including those 
that lie in the low-density regime where a lower-energy 
UHF solution exists (i? > 6 for Z? = 3). 1211 



III. EXPLICITLY CORRELATED CALCULATIONS 

We write the correlated wave function of ZJ-ballium as 

(jj n uj 

* = E E E ^«'™ + ^12) i'nlm, (20) 



n=0 1=0 m=0 



where P12 is the permutation operator between electron 
1 and 2, which ensures the correct symmetry for the 
ground state, and the basis functions are 



^„i™ = (i-a;')(l-y")x^V'^'", 



(21) 
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TABLE I. Hartree-Fock and exact energies of 3-ballium for R — 1, 5 and 20. 



Thompson and Alavi 
Present work using f9| 



Basis set size 

7 
7 
10 



Hartree-Fock energy 



R=l 
11 Ml 747 645 
11.641 747 631 859 
11.641 747 631 855 851 828 



R = 5 
0.739 761 807 
0.739 761 794 626 
0.739 761 794 625 138 



7? = 20 
0.105 378 511 
0.105 378 488 
0.105 378 488 024 



Exact energy 



Thompson and Alavi |^ 

Jung and Alvarellos 
Present work using (20| 



^ Reference |.21j 
^ References [161 and 1311 
Reference ll7l 



R=l 

210 11.591 380 285 

Extrap. 11.590 81(4) 

6296 11.590 906 

196 11.590 838 69 

726 11.590 838 689 02 



R = 5 
0.701 706 934 
0.701 606 (2) 

0.701 613 820 
0.701 613 820 002 



i? = 20 
0.086 577 117 
0.086 577 0(0) 

0.086 576 568 4 
0.086 576 568 358 529 



the scaled coordinates are 

y = 



R' 



1j2 

R' 



ri2 
R ' 



(22) 



and n, I and m are non-negative integers. Such functions 
ensure that 5* is smooth at ri = and r2 = 0, i.e. 



5* 

dri 



ri=0 



94- 

dr2 



0, 



(23) 



r2=0 



and that is cusped at the boundary and satisfies (|4|. 
The total number of basis functions in (21 ) is 



N = 



(a; + l)2(w + 2) 



(24) 



The ground-state energy is the lowest eigenvalue of 

g-i/2(-rj,^U)S-i/2, (25) 

where S, T and U are the overlap, kinetic and repulsion 
matrices, respectively.'^ 

Although our (x, y, z) coordin ates are equivalent to the 
(s, t, u) coordinates of Hylleraas,'^'^ ours lead to simpler 
closed-form expressions for the required integrals. All 
the required matrix elements can be found in closed form 
using the general formula in Appendix \K\ 

Figure [2] shows how the accuracy of the exact energy 
of 3-ballium improves as the number N of terms in the 
expansion increases and Table [T] reports numerical values 
of the exact energy for various i? = 1, 5 and 20. Our 
explicitly correlated results are comp ared with the CI 
energies of Thompson and AlavP^Hll and of Jung and 
Alvarellos. 23 Convergence for other values of D is similar. 

Explicitly correlated calculations converge much faster 
than CI calculations because the former include terms 
(21 1 with m — 1 satisfying the Kato cusp condition.l^^ 
For example, for the unit ball, Thompson and 
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FIG. 2. Accuracy A of the exact energy of 3-balhum with 
respect to the basis set size for various R. TA results taken 
from Ref. 1161 JA from Ref. ilTjand LG from the present study. 



obtained E = 11.591 380 285 using 210 basis functions, 
Jung and Alvarello^I^l subsequently found E = 11.590 
906 using 6296 functions, but we obtain = 11.590 838 
689 using only 196 functions (w = 6). Our energy is 
consistent with the extrapolat ecP^ estimate E = 11.590 
81(4) of Thompson and Alavi.EsEl Likewise, using 6296 
basis functions, Jung and AlvarelloJ^^ found E = 22.033 
71 for the first excited S state and, using 196 explicitly 
correlated functions, we obtain E = 22.033 562 4 E^. 

Figure [2] reveals that, for i? = 1 and R = 5, the rate 
of convergence of the Hylleraas basis set is very similar. 
When R exceeds the Wigner-Seitz critical value {rg « 6 
for D = 3), a Wigner molecule is formed, characterized 
by a m inimu m of the electron density at the center of 
the box.'iSP^' Although the CI calculations of Thompson 
and Alavi show only small variations of the rate of con- 



4 



vergence, the Wigner molecule formation dramatically 
modifies the energy convergence of our explicitly corre- 
lated calculations. For R = 20 and short expansions 
{N < 50) , the CI and Hylleraas calculations lead to simi- 
lar results but, for larger basis sets, the Hylleraas scheme 
is superior and the accuracy rapidly improves. 

Correlation energies for D = 2,3, ... ,7 and R = 
1, 5 and 20 are given in Table |llj For fixed D and in- 
creasing R, Ec decreases. For fixed R and increasing D, 
although both the exact and HF energies increase, Ec 
decreases .'^"'^ 



A. Hartree-Fock energy 

(2) 

Values of E^-p have been determined using the gener- 
alization of the Byers-Brown-Hirschfelder equation^^ 

4f--/' oZ^'i^'f ,2 dr, (34) 
W{r) = 2 / [J^, (x) - Mxf dx, (35) 
where J^,, (x) is given by ([8| . 



IV. LIMITING CORRELATION ENERGY 

Following Hylleraas perturbation theory,'^ we expand 
both the exact and HF energies as series in R, yielding 



E=^ + ^+E(-^+0{R), 
i^HF = + + 4f + 0{R). 
The limiting correlation energy is then given by 



lim - - 41 



The one-electron Hamiltonian for Z3-ballium is 







r rf2 D-l d 

1 



V{r), 



r dr 

and the associated zeroth-order wave function is 



*o(ri,r2) = ^/'o(ri)V'o(?'2), 



where 



-J2 Jd/2-i{i^t) 



(26) 
(27) 

(28) 
(29) 
(30) 
(31) 



In (31), K = jD/2-1,1 and j^_fc is the fc-th zero of the 
Bessel function of the first kind j£)/2-i-'^ The E^^'^ and 
values are easily obtained from the relations 



E 



(0) - ^2 



(32) 



and are reported in Table |ll| For odd D, E^^^ can be 
found in closed form. For example, for D = Z, 



1 



Si(27r) Si(47r) 



27r 47r 
where Si is the sine integral function.!^ 



(33) 



B. Exact energy 

The second-order energy which minimizes the 

Hylleraas functional,'^ is given by 



where 



A = T-£;(")s, 



e'^^^s - u 



(36) 

(37) 
(38) 



The matrices S, T and U have been defined in Sec. |III| 
The vector C contains the coefficients of the zeroth-order 
wave function (30) expanded in the basis set (21). The 



basis set has been enlarged by progressively increasing 
the maximum value of uj. 



C. Correlation energy 

The exact and HF second-order energies, as well as 
the limiting correlation energy E^°\ of i:>-bamum are 
reported in Table |llj The latter is compared with previ- 
ously reported result^^II for related two-electron systems 
(helium, spherium and hookium). 

The first observation is the tendency of the lim- 
iting correlation energies to decrease with increasing 
dimensionality.'^ As D increases, all of the energies de- 
crease dramatically and the correlation energies fall by al- 
most two orders of magnitude between D — 2 and D = 7. 
This point has been already discussed and explained in 
previous works. 

We have used the method developed by Herschbach 
and collaborator J® ^^^^^ to find that the large-_D expan- 
sion of i?c°^ in Z3-ballium is 



1 



53_, 
128^ 



(39) 



where S = 1/(D — 1).!^ This supports our recent con- 
jecture ([T]) that the leading term —(5^/8 is universal and 
independent of the radial external potential V{r). We 
note that the coefficient of 5^ is larger than those in 
the other two-electron systems and this explains why the 
limiting correlation energy of Z)-ballium is always larger 
than those in helium, spherium and hookium. 
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TABLE II. Exact, HF and correlation energies for various finite R. Zeroth-, first-, second-order energies and limiting correlation 
energies. 



D 


2 


3 4 5 6 


7 


E 

E, 


8.I0493I 
8.326496 
-0.221565 


Exact, HF and correlation energies of ballium for R = 1 
11.590839 16.151742 21.519813 27.612654 
11.641748 16.172654 21.530902 27.619443 
-0.050909 -0.020912 -0.011089 -0.006789 


34.391191 
34.395746 
-0.004555 


E 

i?HF 


0.586796 
0.711077 
-0.124281 


Exact, HF and correlation energies of ballium for R = 5 
0.701614 0.863437 1.063334 1.296170 
0.739762 0.880997 1.073192 1.302407 
-0.038148 -0.017560 -0.009858 -0.006237 


1.559045 
1.563317 
-0.004272 


E 

E, 


0.078628 
0.123044 
-0.044416 


Exact, HF and correlation energies of ballium for R — 20 
0.086577 0.096381 0.108129 0.121815 
0.105378 0.107060 0.114985 0.126552 
-0.018801 -0.010679 -0.006856 -0.004737 


0.137388 
0.140835 
-0.003447 


ew 


5.783186 
2.596157 


Zeroth- and first-order energies of ballium, from Eq. ( 32 1 
9.869604 14.681971 20.190729 26.374616 
1.786073 1.496754 1.343463 1.246845 


33.217462 
1.179626 


„(2) 


-0.324120 
-0.057959 


Second-order energies of ballium, from Eqs. (|34[) and (|36|l 
-0.069618 -0.028107 -0.014770 -0.008977 
-0.014442 -0.006194 -0.003333 -0.002037 


-0.005983 
-0.001352 


Helium 
Spherium 
Hookium 
Ballium 


-0.220133 
-0.227411 
-0.239641 
-0.266161 


Limiting correlation energies Sc"' , from Ref. [23 and Eq. 1 28 1 
-0.046663 -0.018933 -0.010057 -0.006188 
-0.047637 -0.019181 -0.010139 -0.006220 
-0.049703 -0.019860 -0.010439 -0.006376 
-0.055176 -0.021913 -0.011437 -0.006940 


-0.004176 
-0.004189 
-0.004280 
-0.004631 



CONCLUSION 



to non-spherical external potentials. 



In this article, we have reported accurate results for 
the exact and HF ground-state energies of two electrons 
of opposite spin confined within a ball of radius i? in a 
ZJ-diniensional space. We call this system D-ballium. 

Our results, focussing mainly on the high-density 
regime (small-i?) and the limit where R = 0, e xtend and 
complete earlier studies on 3-ballium.l^^'^^El] xhe exact 
results have been obtained using a Hylleraas-type expan- 
sion, while the HF limit has been reached using a simple 
polynomial basis set. 

We have also shown that, in the large-dimension limit, 
the limiting correlation energy behaves as —(5^/8 — C6^, 
in agreement with our recent conjecture,'^ and is larger 
than the limiting correlation energy in other two-electron 
systems. A rigorous proof of the conjecture will be re- 
ported elsewhere,'^ where we will show that this result is 
related to the Kato cusp factor. Using continuity argu- 
ments, it seems clear that the conjecture does not apply 
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Appendix A: Hylleraas-type integrals 



The integrals needed to compute the different matrix 
elements are of the form 



(Al) 



6 



with the volume element 



0-3 

dT = xyzJ 2 dxdydz, 



(A2) 

J = {x + y + z){x -y + z){x + y- z){x -y- z), 

(A3) 

and domain of integration 

//•I i-l rx+y 

dT= dx dy dz. (A4) 
Jo Jo J\x-y\ 



One eventually finds 



and 



r( 


Z5-1 
2 


) ]^n+l+m+2D 


r 


(f) 


n + l + m + 2D 



{iT+in, (A5) 



p /a+Z) _b fc+D-2 . a+D+2 D.'i\ 
3 2 1,2' 2' 2 ' 2 '2' / 



(A6) 
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